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ABSTRACT 

We report medium resolution VLT ISAAC K-band spectroscopy of the nuclei 
of seven ultraluminous infrared galaxies. After accounting for stellar absorption 
features, we have detected several molecular hydrogen (II2) v =1-0, 2-1, and 3-2 
vibrational emission lines, as well as the HI Br7 and Hel 2^P — 2^S recombination 
lines. The relative H2 line intensities show little variation between the objects, 
suggesting that the H2 excitation mechanisms in the nuclei are similar in all the 
objects. The 1-0 emissions appear thermalised at temperatures T ~ 1000 K. 
However, the 2-1 and 3-2 emissions show evidence of being radiatively excited by 
far-ultraviolet (FUV) photons, suggesting that the H2 excitation in the ULIRGs 
may arise in dense photon dominated regions (PDRs). We show that the line 
ratios in the nuclei are consistent with PDRs with cloud densities between 10^ 
to lO^cm"^, exposed to far ultraviolet (FUV) radiation fields at least 10^ times 
more intense than the ambient FUV intensity in the local interstellar medium. We 
have constructed starburst models for the ULIRGs based on their H2 properties, 
as well as on the intensities of the recombination lines. Our models provide a 
consistent picture of young 1-5 Myr star clusters surrounded by relatively dense 
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PDRs which are irradiated by intense FUV fluxes. Comparison to the inner few- 
hundred parsecs of the Milky Way indicates that the star formation efficiency in 
ULIRGs is 10-100 times higher than in the Galactic Center. 

Subject headings: molecular processes — galaxies: ISM — galaxies: starburst — 
infrared: galaxies — galaxies: nuclei — line: formation 



1. Introduction 

In star forming galaxies, molecular hydrogen (H2) vibrational emission lines include 
some of the brightest lines in the near-infrared (1.9-2.5 /xm) K-band, and have therefore been 
candidates for detailed studies. The H2 lines may be used to probe the physical conditions of 
molecular clouds near the star forming regions. For example, the relative strengths of lines 
in the v — 1 — and v — 2 — 1 bands allow radiative fluorescent excitation in low-density 
photodissociation regions (PDRs) (Black and van Dishoeck 1987; Sternberg 1988; Sternberg 
and Dalgarno 1989; Burton et al. 1990b; Drainc and Bertoldi 1996) to be distinguished from 
collisional excitation in shocks or dense PDRs (Sternberg and Dalgarno 1989; Burton ct al. 
1990a; Chang & Martin 1991; Draine and Bertoldi 1996; Lim et al. 2002). Determining the 
H2 excitation mechanism in a given source generally requires the observation of many, often 
weak lines, particularly since different processes (e.g. shocks vs. FUV heating) can give rise to 
similar intense and thermalised 1-0 emissions. This requirement has resulted in a great deal 
of uncertainty in the interpretation of H2 emissions in extragalactic objects, because signal 
to noise limitations, the presence of bright and uneven stellar continua, and - particularly 
in earlier work - low spectral resolution, have meant that often only the few strongest lines 
could be detected. 

Recent results from Gilbert et al. (2000) and Davies et al. (2000) show that these 
difficulties can in fact be overcome, and that detailed studies of the molecular hydrogen 

in extragalactic objects is possible as long as there is sufficient signal to noise to detect 
some of the weaker lines. Gilbert ct al. (2000) observed individual star clusters in the 
Antennae (NGC4038/9), which lies at a distance of only 19Mpc. Near the embedded star- 
cluster located in the mid-lR peak they detected 22 H2 transitions, including some from very 
high vibrational levels, extended over an area twice the size of the continuum and nebular 
regions. Gilbert et al. concluded that the emission fines are produced by pure radiative FUV 



^Based on observations at the European Southern Observatory VLT (67.B-0417). 
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fluorescence, most likely in PDRs illuminated by the young star clusters. In work on the 
two nuclei in Mkn266, which at 115 Mpc is too distant to permit studies of individual star 
clusters, Davies et al. (2000) used stellar templates to subtract the continuum, allowing the 
fluxes of the weaker lines to be measured. They found that in neither case could single shock 
or PDR models fit the H2 line ratios. They concluded that ~ 70% of the prominent 2.12 fim 
1-0 S(l) line originates in non-dissociative shocks; and, based on morphology and energetics, 
that the remaining line luminosity arises in X-ray irradiated gas, FUV-pumped gas, and in 
dissociative shocks. 

In this paper wc study H2 excitation in ultraluminous infrared galaxies (ULIRGs). Such 
objects often have 1-0 S(l) to Br7 line ratios greater than unity (Goldader ct al. 1995; Murphy 
et al. 2001). In star forming regions where the 1-0 emission is produced by FUV-pumped 
fluorescence the 1-0S(1)/Br7 intensity ratio is < 1 (e.g. Puxley et al. 1988; Doyon et al. 
1994; Davies et al. 1998), suggesting that in ULIRGs the 1-0 emission is coUisionally excited 
in hot gas leading to higher H2 line surface brightnesses. Indeed, because the majority of 
ULIRGs show signs of elevated star formation rates, as well as interactions and supcrwinds 
(Sanders & Mirabel 1996), one might expect shocks or dense PDRs to play important roles 
in the excitation and emission processes. We have obtained spectra deep enough to detect 
V = 2-1 and 3-2 transitions, with the goal of constraining the H2 excitation mechanisms in 
the ULIRGs. 



2. Observations and Data Reduction 

Our data were obtained as part of the service observing programme at the Very Large 
Telescope between June and September 2001. ISAAC was used in Medium Resolution mode 
in the K-band, providing a nominal resolution of R=2600 with a 1" slit. Since the wavelength 
coverage in this mode is only 0.122 /xm, two settings were used with an overlap of 5-10%. 
The total integration time was 20 minutes per setting. Standard calibrations were performed, 
including atmospheric standard stars (type B or GOV), arcs, flatflelds, and dark frames. The 
data were reduced using PG-IRAF 2.11.3 using standard techniques. Wavelength calibration 
was achieved in two steps: the arc frames were used to apply a transformation which re- 
moved curvature and provided an initial calibration; then the sky lines on the object frames 
themselves were used to remove any residual tilt and flx the flnal calibration. 

Spectra were extracted over 5 pixels (0.75"). Although in a number of cases the line 
emission was extended (in two cases very much so) beyond this, a short aperture was used 
for two reasons: (1) the slit position angles were assigned rather arbitrarily, and so to avoid 
biassing the analysis, we restricted the spectra to the nuclear region; and (2) the signal-to- 
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noise was invariably reduced if the spectra were extracted over a longer aperture, most likely 
as a result of the bright continuum. The two segments of the spectrum for each object were 
combined after applying a scaling measured from the overlapping region. The same scaling 
was then also applied to the reduced 2D spectra. 

An approximate flux calibration was derived from the standard stars, the K-band mag- 
nitudes of which were estimated from the spectral type and V-band magnitude; there was a 
very good correspondance (to better than 0.1 mag) to the K-band magnitudes that could be 
measured directly (e.g. from the 2MASS catalogue). The K-band magnitudes are typically 
up to ~1 mag fainter than values published elsewhere. This is a result of the smaller aperture 
used here (e.g. Goldader et al. 1997 used 3x3 or 3x9 arcsec apertures). Similarly the 1-0 S(l) 
fluxes are typically a factor ~ 2 fainter than those of Goldader et al. (1997). Notably, the 
flux is 3 times fainter in IRAS 16164-0746 and MCG -03-04-014, but position- velocity 
diagrams for these objects show the line emission is very extended; on the other hand, the 
flux is only 50% less for IRAS 14378-3651. 

Acquisition images with integration times of a few seconds were taken before each spec- 
trum. These are shown in Fig. 1 together with the orientation of the sht across the nucleus 
of each object. The position angle of the sht was set to he along the major axis of the object 
as determined from 2MASS K-band or red Digitised Sky Survey images. Where there was 
no clear axis, the slit was set east-west. 

3. Continuum Fitting and Line Extraction 

The weakest v — 3-2 lines that we have detected have strengths that are only a few 
percent of the 1-0 S(l) hne. These are difficult to detect because they are similar in scale 
to the absorption features present in the stellar continua dominated by late type giant or 
supergiant stars. In particular the 3-2 S(5) line at 2.066 ;um coincides with a bump on the 
continuum at 2.067 /xm; and the 3-2 S(3) line at 2.201 fim sits on the edge of the Na absorption 
at 2.207/im. To overcome these difficulties wc have fitted and subtracted a continuum 
constructed from the stellar templates of Wallace and Hinkle (1997). This library includes 
over 100 K-band spectra of stars spanning spectral types 0-M and luminosity classes I- 
V, at a resolution of R~3000. We tried a variety of combinations of template stars, and 
concluded that while there was very little difference between the fits found using giants or 
supergiants, the dwarfs produced a noticeably poorer match to the continuum. Finally we 
used 7 templates of supergiants ranging in spectral type from F5 to M5. Extending the range 
to earlier types or increasing the number of templates used had no noticeable effect on the 
final result. During the fitting process, the relative scalings of these were varied to minimise 
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the difference between their sum and the spectrum, excluding any regions close to line 
emission. The resulting combination of templates was then reddened (i.e. effectively allowing 
adjustments to the overall slope of the continuum) using the Howarth (1983) extinction curve, 
and finally convolved with a Gaussian to take account of the stellar velocity dispersion. The 
extinction and Gaussian FWHM were allowed to vary during the minimisation. In all cases 
except for IRAS 16164—0746, the extinction derived in this way was small {Ay < 1), a result 
that is not surprising if the stars and dust are mixed. The results of the fitting procedure are 
shown in Fig. 2 as the thick lines drawn over the spectra: the fits are successful and in most 
cases match the observed continuum features well. This confirms that the scaling between 
the two segments of each spectrum is correct, and that the weak line fiuxes are likely to be 
accurate. 

Line fiuxes were measured from the spectra after subtraction of the stellar continuum, 
and are given in Table 4. We have detected the S(2), S(l) and S(0) fines in the v — 1-0 
band, the S(3), S(2), and S(l) lines in the 2-1 band, and the S(3) line in the 3-2 band. The 
transition wavelengths are listed in Table 4. We have also set upper limits for the 3-2 S(5) 
line intensities. In addition, in each galaxy we detected the HI 2.17 fan Br7 and the He I 
2^P — 2^S 2.06 fim recombination lines. The uncertainties quoted are measured directly 
from the spectra as the rms of the residual and hence are rather higher than the errors that 
would be estimated purely from photon noise, including as they do the effects of improper 
line and continuum fitting. Based on this, we have reached a mean signal-to-noise of 25a for 
the 1-0 S(l) line. 

It is apparent from Table 4 that although there is considerable variation in the ratios of 
He I and Br7 to 1-0 S(l), the H2 line ratios are in all cases very similar. Given the range of 
values which the ratios could take amongst these objects, the standard deviation in the ratio 
of any particular H2 line to 1-0 S(l) of only 0.02-0.03 is small. We have therefore derived a 
mean set of ratios, which we also include in Table 4. 



4. H2 Excitation Diagrams and PDR Models 

Given the H2 line fiuxes / (W m"^) listed in Table 4, we have calculated the impfied 
molecular column densities, N^j = ATrf/AQ,, in the upper rotational-vibrational (vj) levels 
of the observed transitions, where A is the radiative rate and Q — 0.75"^ is our aperture 
size. We assume that the quadrupole transitions arc optically thin and we use the radiative 
A values as given by Wolniewicz et al. (1998). In Figs. 3 and 4 we show plots of log Nyj / g^j 
(where g^j is the statistical weight) vs. the energy (in K) of each observed vj level. We have 
normalized the population distributions relative to iVi a/gii 3 as inferred from the 1-0 S(l) 
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line. 

For thermalised populations at fixed gas temperature the logN/g points should he on 
a straight line in the excitation diagrams. Clearly, such single component models do not fit 
the data. While the v — 1 levels appear to be thermalized at T ~ 1300 K, inclusion of the 
V — 2 levels raises the best-fitting temperature to T ~ 2400 K. Even at this temperature, 
however, the v = 3 levels remain underpopulated, and the excitation temperature inferred 
from the relative populations in the v = 2 and v = 3 levels exceeds 5000 K. At such gas 
temperatures the molecules would be rapidly dissociated, suggesting that a non-thermal 
excitation mechanism is responsible for the excitation of the v = 2 and v = 3 levels. 

Collisional excitation in shocks could possibly account for the 1-0 emissions. A line 
intensity of 10"'^ W m^^ sr^^ is equivalent to a flux of 2 x 10"^^ W m~^ for our 0.75"^ aperture, 
consistent with the fluxes given in Table 2. It is only in models of non-dissociative C-shocks 
with velocities > 25kms~^ and preshock densities > lO^'cm"^ that the predicted 1-0 S(l) 
intensity reaches such values (Burton et al. 1990a; Timmermann et al. 1998; Wilgenbus et 
al. 2000). In C-shocks at speeds lower than this, and in all realistic J-shocks, the 1-0 S(l) 
intensity is significantly lower. In C-shocks at or above this speed (and for gas densities 
> 100 cm~^) the post shock temperatures are of order 1000 K. Such models would have 
difficulty in accounting for the relatively large columns in the u = 2 and i/ = 3 levels. 
Furthermore, for shocked emission the ortho-to-para ratio in vibrationally excited levels 
should attain an equilibrium value of 3 (Timmermann et al. 1998; Wilgenbus et al. 2000). 

In contrast to shocks, the observed H2 excitation could be produced in dense PDRs. 
In such clouds, FUV irradiation can heat the gas to temperatures ~ 1000 K, leading to 
collisionally excited and thermalized 1-0 line emissions, while direct FUV-pumping main- 
tains non-thermal fluorescent populations in higher-lying vibrational levels (Sternberg and 
Dalgarno 1989; Burton et al. 1990b; Draine and Bertoldi 1996). In fact, a high vibrational 
excitation temperature, such as the 5000 K inferred from the relative v = 3 and v = 2 
populations, is a basic characteristic of FUV-pumped fluorescent emission (Black and van 
Dishoeck 1987; Sternberg 1988). Furthermore, an additional signature of fiuorescent excita- 
tion appears to be present in our data set. For those objects in which the v — 2 — 1 S(l), 
S(2), and S(3) lines were detected, the resulting N/g values for the j = 3, 4, and 5 levels do 
not appear thermalised. Instead, the j — 4 para-level is shifted above the adjacent j = 3 
and J = 5 ortho-levels. Thus, the ortho-to-para ratio for the v = 2 molecules appears to 
be out of equilibrium and suppressed below the ortho-to-para ratio of 3 obtained in local 
thermodynamic equilibrium (LTE) in warm gas. This behavior is precisely what is expected 
for excitation by FUV-pumping, in which the greater optical depths in the ortho UV ab- 
sorption lines supresses the populations in the vibrationally pumped ortho-H2 relative to 
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para-H2 (Sternberg and Neufeld 1999). 

We have constructed a series of representative PDR models to determine PDR conditions 
and parameters consistent with the observed H2 hne ratios. For this purpose we have used the 
code described in Sternberg and Dalgarno (1989, 1995), and Sternberg and Neufeld (1999). 
Briefly, our models consist of static, plane-parallel, semi-infinite clouds that are exposed to 
isotropic FUV radiation fields. At each cloud depth we compute the equilibrium atomic to 
molecular hydrogen density ratio, n{H)/n{H2), and we solve for the steady-state population 
densities in the rotational and vibrational H2 levels in the ground electronic state. In solving 
for the vj populations we include the effects of FUV-pumping via the Lyman and Werner 
bands, coUisional processes with H"*", H, and H2, and quadrupole radiative decays. 

We present a series of five representative PDR models in which we consider a range 
of illuminating FUV field intensities X) relative to the FUV field in the local interstellar 
medium (2.1 x 10^^ photons m~^, Draine 1978), and total hydrogen particle densities 
nu = n(H) + 2n(H2) in the clouds. In all of the models except model 1 (for which we keep 
the temperature constant at T = 100 K), we assume that the gas temperature varies with 
cloud depth as 

T 

rp max /-|^\ 

~ 1 + 39 X (2n(H2)/nH)5 ^ ' 

where Tmax ~ 10^ K is the temperature at the cloud edge. For such thermal profiles the gas 
temperature equals T^^ax in the outer atomic zone of the PDR, and then declines significantly 
to ~ 25 K as the gas becomes molecular. This behavior is consistent with theoretical 
expectations for dense PDRs, and is also empirically based (e.g. the Galactic star forming 
cloud S140, as analyzed by Timmermann et al. (1996) and Sternberg and Neufeld (1999)). 

We summarize the parameters for our five models in Table 3. The level populations from 
the models are given in Table 4 (below the observed values), and are plotted in Figure 3. 
Also shown in the figure are the populations derived from the mean observed line ratios. We 
now give a short summary of the models and how well they describe the mean observed line 
ratios: 

model 1 Low density, nn = lO^cm"^, and cool (isothermal) T — 100 K gas, excited by a 
relatively weak x = 10^ FUV field. This results in low surface brightness H2 emission, 
with /i_os(i) = 1-3 X 10~^ W m~^ sr~^. The population ratios, including the v = 1 levels 
are characteristic of pure radiative fiuorescent excitation, unaffected by vibrational 
collisional excitation or deexcitation. This model does not produce a good fit to the 
data. 

model 2 The gas density and UV field are both increased by an order of magnitude to 
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riH — lO^cm"^ and x — 10^; ^ thermal profile as given by Eq. 1 is adopted 
with Tmax = 10^ K. This results in much higher line surface brightnesses, especially 
for the 1-0 transitions, with /i_os(i) = 5.2 x 10"*^ W m^^ sr~^. The v = 1 levels are 
thermalised by collisions, while the f = 2 & 3 levels are excited by FUV-pumping 
and therefore display the ortho-to-para shifts characteristic of fluorescent excitation 
(Sternberg and Neufeld 1999). This model matches the data rather well, showing that 
the H2 populations are consistent with those expected for dense PDRs with hot outer 
boundaries. 

model 3 Similar to model 2, with again nn = 10"^ cm^^ and x = 10^, but with T^ax increased 
to 2000 K. This increases the surface brightnesses in the 1-0 lines still further, with 
-^i-os(i) = 2.0 X 10~^ W m~^ sr~^. The higher v levels are now also thermalised and 
therefore have much lower populations relative to the v — 1 level. The result is very 
similar to a simple purely "thermal" model in which the populations are assumed to 
be in LTE at 2000 K, and very much under-predicts the v — 2 3 populations. 

model 4 Also similar to model 2, with nn = lO^cm"^ and T^ax = 10'^ K, but with the 
UV field increased by a factor of 100 to x = 10^. The relative populations and hne 
intensities are hardly affected compared to model 2, and there is only a slight increase 
in the populations of the higher v levels. 

model 5 As in model 4, having x = 10^ ^iid T^ax = 10'^ K, but with the gas density 
increased by a factor of 100 to nn = 10^ cm^^. The higher density leads to more 
effective thermalisation in the higher v levels, so that these levels are under predicted 
by the model, although not as severely as in model 3. 

It is clear that models 2 and 4 provide the best match to the mean line ratios, and that 
with the current data they are effectively indcstinguishablc. The 1-0 S(l) line intensities of 
the 2 models are also similar, and consistent with those observed in the ULIRGs. 

In our models we have not explicitly computed the heating and cooling balance that 
would be expected to yield hot gas in the outer photodissociated atomic zones. In PDR 
theory the 6-13.6 eV photons are available for heating via photoelectric emission from dust 
grains, and collisional deexcitation of FUV-pumpcd H2. For our model 2, the total energy 
flux in this band is 2.6 x lO^'^Wm"^; while the emergent flux in the II2 1-0 S(l) line 
is 5.2 X 10~®Wm~^, or 2 x 10~^ of the available FUV energy, a fraction consistent with 
theoretical computations (Sternberg and Dalgarno 1989; Burton et al. 1990b). Thus it is 
possible that the H2 emission - including the thermalised v — 1 — emissions - is produced 
entirely in PDRs. However, for the 1-0 lines we cannot rule out contributions from additional 
sources such as shock waves. 
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Having looked at which models might match the mean ratios, we now briefly consider 
the individual objects. 

Prom surface brightness considerations alone, model 1 can be ruled out for all objects. 
Integrated over the 0.75"^ aperture in which we have measured the line fluxes, this model 
would predict a maximum 1-0 S(l) flux of 2.3 x 10'"^'^ Wm~^, nearly 2 orders of magnitude 
less than the range observed in our sample. We would also argue against model 3 because 
it requires that all the levels are effectively thermalised at the same temperature. We now 
consider the remainder of the models for each object. For IRAS 19458+0944, models 2 & 4 
provide good fits to the data, and while model 5 also fits it is more to the extreme of the 
errorbars. These three models also match the data for MCG — 03— 12— 002 well, although 
the V — 1 levels indicate a temperature at the edge of the clouds somewhat warmer than 
1000 K. For MCG -03-04-014 models 2 and 4 also match the data weU, but here the ex- 
tra data point for v = 2 levels suggests that model 5 is less hkely. IRAS 01364— 1042 and 
IRAS 16164—0746 present some of the best data. For these objects model 5 van be con- 
fidentally dismissed in favour of models 2 and 4. For the former, the v = 1 populations 
suggest, as for MCG —03—12—002, that the temperature at the edge of the clouds is slightly 
warmer than 1000 K. The cases for IRAS 14378-3651 and IRAS 20414-1651 are rather un- 
certain and models 2, 4, and 5 all match the data well; although in the latter object, there 
is a tendency towards the higher density model since the v — 2 distributions appear more 
thermalised. 

Our conclusion from the line intensities and ratios is that all the objects can be charac- 
terised by simple single-component PDRs illuminated by FUV fields with x ^ 10'^, and hy- 
drogen densities in the range nn = 10^-10^ cm~^ (possibly in some cases as high as 10^ cm~^), 
and a temperature at the outer edge of the clouds T — 1000-1500 K. Multi-component mod- 
els (Burton et al. 1990b) could also be considered, in which, for example, an inner dense 
region is responsible for the thermalised 1-0 emissions, with contributions to the fluorescent 
component from more extended low density gas (see also our discussion of Orion in §4.2). 
However, given the limited spatial information in our data set we do not consider such models 
here. 

Observations of additional FDR diagnostic lines such as [O II] 63 /xm and 145 fim, [C II] 
158 fim, or high-J CO transitions (Tielens and Hollenbach 1985; Sternberg and Dalgarno 
1989; Koester et al. 1994) would provide valuable additional constraints on the PDR condi- 
tions. However, no such data is yet available for our sources. 
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4.1. Starburst Models 

If the observed molecular hydrogen emissions are produced in PDRs then it is of interest 
to consider whether the FUV sources - most likely hot OB stars - can also provide the Lyman 
continuum (Lyc; hu > 13.6 eV) EUV photons responsible for the Br7 recombination line 
emission. As we will now show, the relative photon production rates in the FUV and EUV 
bands that are implied by the H2 and Br7 line fluxes, are consistent with young star clusters 
containing populations of massive hot stars. 

In Table 6 we present the results of starburst model computations using our starburst 
code STARS (e.g. Sternberg 1998; Thornley et al. 2000) for the Lyc and FUV-band photon 
rates Qhyc and produced in evolving star clusters. We deflne the FUV-band as spanning 
912 to 1130 A, since photons in this wavelength range can be absorbed in the H2 Lyman 
and Werner absorption line systems, leading to FUV-pumped fluorescent H2 emission lines. 
In Table 6 we also list the helium ionizing (He I; hu > 24 eV) photon production rates, 
and the associated K-band himinosities, Lx ^, and bolometric luminosities Lboi of the model 
clusters. We consider a range of chister ages, for both instantaneous (10^ ^0) bursts and 
continuous {SFR = 1 Mq yr~^) star formation, assuming Salpeter initial mass functions 
(IMFs) with upper-mass limits of M^p — 30 and 120 Mq. Our starburst models incorporate 
hot-star "wind-atmospheres" in the computation of the cluster luminosities and spectral 
energy distributions (Sternberg et al. 2003). We assume solar metaUicity for the evolutionary 
tracks and stellar atmospheres. 

Table 6 shows that the ratio (5H2/QLyc varies with cluster age and also depends on M^p. 
For instantaneous bursts, Qn2/Q'Lyc increases rapidly with age as the hottest Lyc producing 
O stars disappear from the system, whereas FUV photons continue to be produced by the 
B-stars. For continuous star formation Qn2/Q'Lyc increases slightly with time and reaches 
flxed values as the relative numbers of O and B type stars reach equilibrium. At early times 
the photon ratio is smaller for larger M^p due to the efficient production of Lyc photons in 
the most massive stars. 

We determine Qljc for our sources from the Br7 fluxes assuming case B recombina- 
tion at 10*^ K, for which 71.9 Br7 photons are produced for every Lyc photon absorbed in 
photoionizations. Similarly, (5h2 is proportional to the flux in a FUV-pumped fluorescent 
H2 emission line. To estimate (5h2 use the flux in the the 2-1 S(3) line since this is a 
FUV-pumped line that we detected in every galaxy in our sample. The FUV-band to H2 



^we define Lk as the total luminosity in the 1.9-2.5 /zm K-band for a frequency independent K-band flux 
density. If 52.2 is the 2.2 /xm flux density in mJy and .Dmpc is the source distance in Mpc then Lk{Lq) = 
1.17 X lO^Dgj S'2.2. The absolute K-magnitude Mk = -0.33 - 2.51ogLK 
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line photon "conversion efficiency" depends on tfie ratio of tfie incident FUV intensity to 

gas density, x/'^h, a PDR parameter tliat controls the fraction of the incident FUV photons 
absorbed by molecules as opposed to dust grains in the PDRs (Black and van Dishoeck 
1987; Sternberg 1988). The efficiency is largest when xl^n ^ 0.01 cm'^, and decreases for 
larger values of the FUV to density ratio ^. We assume a FUV to 2-1 S(3) conversion factor 
of ?7 = 4.1 X 10"'^, as given by our PDR model 2. For this model 7] is at about half the 
maximum efficiency. It follows that the distance independent ratio 

Q-R2 ^ -^2-18(3) 1 -^2-1 g(3) ^ 3 25 -^2-18(3) 
FBr7 71.9 T] Xsr-y Fbij 

In Table 5 we list the FUV to Lyc photon ratios, that we infer for our ULIRG sample. The 
ratios range from 0.3 to 0.9, with Qr^/Qljc ~ 0.5 a representative value. 

Our starburst models thus show that the H2 and Br7 line emissions may be produced 
by the same populations of star clusters in the ULIRG nuclei. If the hydrogen molecules in 
the PDRs are absorbing the FUV photons at close to maximum efficiency, then the clusters 
must be young, as shown by the photon ratios listed in Table 6. The clusters must then 
also contain massive stars with Mup close to 120 M©. For instantaneous bursts, a ratio 
Q}i2/Qi'yc ~ 0.5 implies an age of ~ 1 Myr. For continuous star formation the limit on 
the age is less severe, and could be as large as ~ 10 Myr, especially if the FUV to H2 line 
conversion efficiency is smaller than we have assumed. It appears then that the same stellar 
radiation sources can plausibly account for both the H2 lines and the recombination line 
emissions. This would not be the case if, for example, the Br7/H2 line intensity ratios were 
much larger than observed. 

The 2.06 //m He I recombination line provides an additional constraint, as this line is 
proportional to Queh the number of He I continuum photons {hu > 24.6 eV) produced per 
unit time. The strong 2.06 //m line in our spectra is a further indication of the presence 
of very hot and massive stars. However, as noted by Shields (1993) converting the helium 
line flux to Quei depends on a variety of largely uncertain parameters for the H II regions, 
including the nebular densities, gas filling factors, and dust content. Here we derive a 
lower limit, Qil^, assuming that all of the He I continuum photons are absorbed in helium 
ionizations, followed by (case-B) recombinations. The effective recombination coefficient for 
the 2.06 urn line (Smits 1991; Osterbrock 1989) ranges from afgg = 4.74 x lO'^^cm^s"^ at 
low densities (10^cm~^), to af^^Q — 6.78 x lO^-'^^cm^s"-'^ at high densities (10^cm~^), where 
coUisional depopulation of the 2^5" level to the 2^5" and 2^P levels (rather than photon 



^The efBciency also depends on the effective FUV grain absorption cross section in the PDRs. We assume 
a cross section equal to 1.9 x 10~^^ cm~^ as appropriate for Galactic dust (Draine 2003). 
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decay to l^S) enhances the He I 2.06 //m intensity. Given a total case-B recombination 
Q;B(HeI) = 2.73 x lO-^^cm^s"^ (Osterbrock 1989), it follows that at least 5.8 2.06 /xm line 
photons are emitted per He I ionizing photon The resulting limits, Qnc'V^Lyc, are listed 
in Table 5. They are consistent with the young ages we have inferred for the starbursts. 

Without any extinction corrections, the Br7 fluxes and source distances imply Lyc emis- 
sion rates ranging from 6x lO^^ g-i (j^ MCG-03-12-002), to 3.4 xlO^^ s"^ (in IRAS 19458+0944) 
For young 1 Myr instantaneous bursts this implies total stellar masses ranging from 6.7 x 10^ 
to 3.8 X 10^ M© for a Salpeter IMF ranging from 0.1 to 120 M©. The implied stellar masses 
increase significantly if such clusters are more evolved (c.f. Thornlcy et al. 2000). For con- 
tinuous star formation, the Br7 fluxes imply star formation rates from 0.7 to 38MQyr^^, 
assuming an age of 5 Myr. It follows from the starburst data in Table 6 that for young 
(IMyr) instantaneous bursts the predicted K-band luminosities of clusters account for at 
most a few percent of the observed K-band luminosities. For continuous star formation, 
the fractions could be a bit larger, up to ~ 20%, but are still small. We conclude that the 
observed K-band continua are not produced in the current star forming episode. 

Equivalently one can consider the Br7 equivalent widths, which should exceed 100 A or 
more for star formation with ages of 10 Myr or less. The directly measured equivalent widths, 
however, lie in the range 5-20 A for these ULIRGs. Low Br7 equivalent widths have also 
been observed by other authors (Goldadcr ct al. 1995; Murphy et al. 2001). The depth of the 
continuum absorption features (in Section 3 we showed that the scale of these features was 
similar to that in late- type giants and supergiants) suggest that it is unlikely that the K-band 
continua arise from hot dust associated either with starbursts or hidden AGN. On the other 
hand, it is plausible that the continua are produced by an older underlying population. The 
only firm conclusion we can draw about such a population is that it contributes little to the 
total FUV and ionising fiux, which means that either it is at a low level or has ceased in the 
recent past. For evolved 10-^° yr stellar populations, M/L^ ~ 20 (Thronson and Greenhouse 
1988; Fig 3 of Tacconi-Garman et al. 1996), implying stellar masses of order 10^° Mq. If the 
near-infrared continuum is dominated by a younger 10^ yr population, consistent with that 
over which star formation occurs in the major merger models of Mihos and Herquist (1996), 
then the mass to light ratio is M/L^ ~ 0.7 (Tacconi-Garman et al. 1996), implying much 
more moderate masses in the range 10^-10^ Mq. In this case, the star formation must have 
ceased already, since it would otherwise make a significant contribution to the observed FUV 
and ionising fiuxes. We note that in some models of Mihos and Herquist (1996) (notably the 
bulgeless mergers), multiple separate epochs of star formation can be triggered during an 



We note that since each He recombination results in about one photon that can ionize H, we have 
assumed that our estimates of Qhyc are unaffected by the presence of He. 
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interaction, occuring at each point of closest approach as the two progenitor galaxies orbit 
each other before finally merging. 

Prom Tables 5 and 6 it can also be seen that the bolometric luminosities based on the 
star formation rates implied by the Br7 fluxes are very much less than those inferred by 
the IRAS fluxes. The total correction factor needed to reconcile them is in the range 20-40. 
The difference could arise from a combination of aperture effects (discussed in Section 2), 
extinction, and perhaps also dust absorption within the nebulae as proposed for ULIRGs by 
Voit (1992) and Bottorff et al. (1998). Mid infrared spectroscopy by Genzel et al. (1998) 
indicates that values of Ay = 5-50 (screen model) or Ay = 50-1000 (mixed model) are 
not uncommon - i.e. for lines in the near infrared one may be seeing as little as 1-10% of 
the total flux. Additionally, dust within the ionised nebulae could absorb the Lyc photons. 
FUV photons would also be absorbed, reducing the FUV intensity at the PDR boundaries. 
Because the dust absorption and aperture effects can apply to the line emitting regions and 
older stellar population differently, it is not straight forward to correct for them and we have 
not attempted to do so. 

Assuming that the observed H2 emission is distributed uniformly over the 0.75"^ aper- 
ture, the area averaged FUV intensities within this region are given by 



where R is the radius probed by our observational aperture. For our ULIRG sample, (x) 
ranges from 700 to 2000, broadly consistent with our PDR model 2. 

The FUV fluxes could be substantially larger if the filling factor of the ionized and 
PDR gas is small. For example, we consider a collection of identical and young star clusters 
each surrounded by its own HII region (which does not overlap with any other), with PDRs 
at the outer boundaries of these HII regions. For such a system, the total Lyc emission 
rate Qhyc = Nc\ns{4:n/3)aBnlrg where A^cius is the number of star clusters, Ug is the electron 
density in the HII regions, the Stromgren radius, and as is the hydrogen recombination 
coefficient. The FUV-band flux at the edges of the PDRs is ((5H2/-^cius)/47rr^, so that the 
FUV flux (in units of the Draine (1978) field) may be written as 



To estimate the mass of a typical star cluster, we consider the sample of 9 starburst 
galaxies studied by Meurer et al. (1995) (see also Maoz et al. 2001). They found that the 
compact star clusters had minimum masses in the range 6 x 10^-2.4 x 10^ M0 assuming a 




(3) 
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Salpeter IMF in the range 5-100 Mq. By extending the lower mass hmit to 0.1 Mq, this range 
increases to 3.3 x 10^-1.3 x 10^ M©. We adopt 10^ as representative of the typical mass. 
For clusters of this mass the Lyc production rate per cluster at 1 Myr is QLyc/-^cius = 9 x 10^^ 
s~^ For Qu^/Qi-yc — 0.5, and ng between 100 and 10^ cm~^ the FUV flux ranges from 2 x 10^ 
to 4 X 10^ 

Wc note that the electron densities in ULIRGs reported by Veilleux et al. (1999) have 
a distribution which peaks at values around 100-200 cm^^ with a long tail to higher values, 
and a few above 1000 cm^'^. However, it is not clear whether these densities, which are 
derived from the [Sll]A6716/A6731 doublet measured in a ~ 2 x 4kpc aperture, do trace the 
density in the inner star forming regions or whether they are dominated by more extended 
emission. In the latter case, the emission could arise in diffuse gas at relatively large radii, 
e.g. in superwinds, and would thus provide only a lower limit on the central density. The 
density profiles for a number of galaxies with superwinds have been measured by Heckman 
et al. (1990). They find that at large radii, greater than a few kpc for the more luminous 
galaxies (which include several ULIRGs), the density is consistent with a r^^ profile and 
has a value of 100 cm~^ at some point typically within r ~ 1-3 kpc. But rig attains a more 
constant value in the range 500-1000 cm~^ at smaller radii. 

In summary we conclude that our Br7 and H2 line fluxes are consistent with nuclear 
star forming clusters containing massive stars, with ages ~ 1 Myr for instantaneous bursts, 
or more evolved 5-10 Myr systems for continuous star formation. The H2 line fluxes, and 
corresponding FUV-band photon production rates, arc consistent with the intense FUV 
fluxes, X = 10^-10^ we adopted in our FDR models 2, 4, and 5. As we have seen, the H2 
line ratios imply that the PDRs are dense, with riu between 10^ and 10^ cm~^, and hot, with 
T = 1000-1500 K in the outer edges. 



4.2. Comparison of PDRs to Other Star Forming Regions 

How do the physical conditions for the PDRs - UV intensity, temperature, density - 
that we have derived in Section 4.1 compare to those in other star forming regions? 

One of the best studied PDRs is that in Orion. Several authors (e.g. Luhman et al. 
1998; Marconi et al. 1998; Usuda et al. 1996) have found that in some regions of Orion 
there is clear evidence for dense PDRs illuminated by intense UV radiation. The integrated 
1-0 S(l) luminosity within a few arcmin of Ori C, which includes all these dense PDRs as 
well as some outflows with shocked emission, is only 4.9 L© (Usuda et al. 1996). However, 
large scale observations of the Orion cloud suggest that the total 1-0 S(l) luminosity is at 
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least 34 and over the entire ~ 3° x 10° CO cloud may be as much as 63 Lq (Luhman et 
al. 1994). That is, the extended emission dominates the total emission by a factor of 7-13. 
Based on the surface brightness (1.4 x 10~^ Wm~^ sr~^) and the strength of the 6-4 Q(l) line 
at 1.601 fim (1-0 S(l)/6-4Q(l)= 2.9) in the extended emission, the authors concluded that 
the excitation mechanism is fluorescence by a weak UV field in low density PDRs, similar 
to our PDR model 1. Hence, although Orion does contain a number of dense intensely 
illuminated PDRs, on a global scale - i.e. if it were observed from far enough away to 
be spatially unresolved - its measured properties would be those of the extended emission. 
Hence, it does not appear to be similar to the PDRs in ULIRGs. By analogy, one would 
might speculate that in ULIRGs there are some PDRs which have much greater (or smaller) 
densities and incident UV fields than the global properties we have been able to measure, 
but that these account for only a small fraction of the total H2 luminosity. 

The extended emission in the inner 400 pc of the Galaxy is different to that in Orion. 
Pak ct al. (1996) found that although there is very intense emission close to SgrA, the 
emission extended on degree scales has a lower surface brightness of ~ 3 x 10"®Wm~^sr~^ 
(corrected for foreground extinction of = 2.5), although this is still an order of magnitude 
higher than that in Orion. Based on this and the uniformity of the emission, the authors 
concluded that it was excited in moderately dense PDRs with moderate UV field - similar 
to our PDR model 2, which has a similar density to but a lower UV intensity than we have 
observed in ULIRGs. Table 3 shows that the intrinsic surface brightness of 1-0 S(l) emission 
is similar in both cases. 

The 400 pc scale over which these measurements were made is similar to that covered 
by our 0.75" aperture at the distance of the ULIRGs. The 1-0 S(l) line luminosities for the 
ULIRGs given in Table 2 has a mean of 2 x 10^ Lq without extinction or aperture corrections. 
Including this correction factor, we find that the line luminosity is 3-4 orders of magnitude 
greater than the 8000 Lq (already corrected for extinction) measured in the Galactic Center 
(Pak et al. 1996). Since we have shown that the surface brightness of the individual PDRs 
is similar both in the Galactic Center and in ULIRGs, the difference in total H2 luminosity 
must be due to a vast increase in the filling factor of the PDRs. 

We argue that this is due not only to an increase in the number of molecular clouds, but 
also an increase in the number of molecular clouds that have become PDRs. The gas mass 
in the inner 500-600 pc of the Galaxy has been determined by a variety of means including 
"^^CO(l-O) and C^^O(l-O) radio measurements, thermal dust emission detected by IRAS and 
COBE, and 0.1-1.0 GeV 7-rays (see Dahmen et al. 1998, and references therein). These all 
yield similar values, for which the authors estimated the weighted mean to be 3 x 10^ Mq. 
Typically, ULIRGs have a strong central concentration of molecular gas on the order of 



-16- 



1O^~^°M0, for example as deduced from observations of the CO(l-O) line by Evans et al. 
(2002) with a 2-4" beam and for a much larger sample, but also with a larger beam of 13- 
22", by Solomon et al. (1997). For the 3 ULlRGs in this paper for which the radio CO 1-0 
luminosities have been measured, the derived gas masses are 1.5-2.0 x 10^° Mq (Mirabel 
et al. 1990). These results indicate that the ratio Li_os(i)/^gas is 10-100 times greater in 
ULIRGs than in the Galactic Center, implying a much higher star formation efficiency for 
the ULIRGs (c.f. Sanders & Mirabel 1996, and references therein). 



5. Spatial Variation of Line Flux 

In Section 4 we have discussed only the line emission in the centre of the ULIRGs. 
Figure 5 shows that in many cases it is extended beyond this, and with a 1-0S(1)/Br7 ratio 
which varies considerably from that in the nucleus. In Figure 6 we show position-velocity 
diagrams for the objects, which show that (where it can be observed) the 1-0 S(l) line has 
the same velocity characteristics as the Br7 line, suggesting that the two lines arise from 
similar regions. In contrast if, for example, the H2 emission arose from material blown out 
or shocked in a superwind (e.g. as for NGC253, Sugai et al. 2003) one might expect to see 
broader or double-peaked profiles similar to those observed in the Ha lines by Heckman et 
al. (1990). However, since the line intensity is very low in the circumnuclear regions, we 
cannot rule out this possibility. 

MCG —03—04—014 appears to be unique among these objects in having the brightest 
Br7 knots offset from the nucleus by 1-2" (0.7-1.4 kpc). This indicates that there is signifi- 
cant recent star formation in the circumnuclear region. At these locations the 1-0S(1)/Br7 
ratio has decreased by a factor of 3 with respect to that in the nucleus, perhaps indicative 
of lower FDR gas densities. 

In many of the other ULIRGs, particularly IRAS 01364-1042 and IRAS 14378-3651, 
there is evidence that the 1-0S(1)/Br7 ratio increases in the circumnuclear region, although 
the surface brightness diminishes rapidly. In IRAS 01364— 1042 the ratio of 1-0 S(l) to the 
continuum peaks at radii of 0.5-1.0" (0.5-1.0 kpc), while the Br7 distribution is very cen- 
tralised. As in MCG —03—04—014, the reduction in 1-0 S(l) surface brightness could be due 
to lower gas densities in the PDRs. IRAS 16164—0746 appears to have very broad extended 
line emission with little change in the 1-0S(1)/Br7 ratio over at least 2" (1.1 kpc). The 
continuum emission is also broad, with a spatial FWHM of 1.4" (0.8 kpc). The image in 
Fig. 1 shows that the slit was aligned close to the position angle of a prominent bar, and 
this may be influencing the observed size of the nucleus. 
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Since the slit position angles for all the objects were assigned rather abitrarily (due to 
the quality of previously available images) and in most cases do not lie along either the major 
or the minor axis of the galaxy, it is difficult to draw strong conclusions about the extended 
emission. However, it is clear that the characteristics of the emission lines (strength and 
ratios) change in the circumnuclear region, and hence the excitation conditions must also 
change. 

6. Conclusions 

We have presented K-band spectra of 7 ULlRGs. Fitting stellar templates to line-free 
parts of the spectra has allowed us to fit and remove the absorption features in order to 
improve our detection limit for faint lines. We have measured line fluxes of Br7 and H2 lines 
from the z/ = 1, 2, and 3 levels, and conclude the following: 

(1) The nuclear H2 line ratios show very little variation between the objects, suggesting 
that the excitation mechanism and conditions are similar in all cases, at least in the nuclei. 
However, the 1-0 S(l) and Br7 line ratio and equivalent widths vary away from the nucleus, 
suggesting that in the circumnuclear region the excitation mechanisms may be different, 
perhaps as a result of a different evolutionary phase of the star formation. 

(2) We find that neither shocks nor models of low density PDRs can satisfactorily fit 
the nuclear H2 line ratios; but FUV-pumped gas in high density PDRs can account for the 
observed ratios and fluxes. Parameter ranges that produce PDR models consistent with the 
observed H2 line ratios are: gas densities riu — 10^~^cm~^, temperatures at the outer edge 
of the clouds T f« 1000 K, and illuminating FUV flelds x = 10^-10^ times more intense than 
the local interstellar held. 

(3) We have constructed a simple model of star formation based on fluorescence of gas 
in PDRs at the edge of HII regions surrounding compact star clusters. This model can 
account for the ratios of the He I, Br7, and fluorescent H2 line intensities, and implies young 
(1-5 Myr) clusters containing massive (120 Mq) stars. Our model implies that the nuclear 
K-band continuum is dominated by a separate older stellar population. 

(4) The PDRs in ULIRGs and those in the inner few hundred parsecs of the Galaxy 
have similar 1-0 S(l) surface brightnesses. Since the 1-0 S(l) emission in both cases is pre- 
dominantly from PDRs, the quantity Li_os(i)/-^gas traces the ratio of the number of PDRs 
to the number of molecular clouds - i.e. the star formation efficiency. This ratio is 10-100 
times larger in ULIRGs than in the Galactic Center. 
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Table 1. Basic Data for the Galaxies Observed 



object 


RA 


dec 


cz 


Sl2 


S25 


Seo 


Sioo 




(J2000) 


(J2000) 


(kms"-') 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


MCG -03-04-014 


01 10 08.9 


-16 51 10 


10040 


0.301 


0.846 


6.48 


10.44 


IRAS 01364-1042 


01 38 52.9 


-10 27 12 


14520 


<0.077 


0.396 


6.16 


6.70 


MCG -03-12-002^ 


04 21 20.0 


-18 48 39 


9652 


0.185 


0.430 


5.75 


8.20 


IRAS 14378-3651 


14 40 59.4 


-37 04 33 


20277 


<0.116 


0.525 


6.19 


6.34 


IRAS 16164-0746 


16 19 11.9 


-07 54 03 


8140 


<0.350 


0.562 


10.20 


13.72 


IRAS 19458+0944 


19 48 15.5 


+09 52 02 


29964 


<0.250 


<0.278 


3.95 


7.11 


IRAS 20414-1651'" 


20 44 17.4 


-16 40 14 


26107 


<0.647 


0.346 


4.36 


5.25 



'^MCG —03—12—002 is a pair with a separation of 17"; we observed only the northern component, 
to which the IRAS data is closer. 

^IRAS 20414—1651 has two nuclei separated by 2.6", of which only one is visible in the K-band. 



Note. — 



IRAS data is from the Faint Source Catalogue and Point Source Catalogue 
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Table 2. Measured Data 



object 




Distance*^ 


K-mag Fsd)^' log^f^'^^ 






(Mpc) 




MCG -03-04-014 


145 


13.4 2.85 6.27 


IRAS 01364- 


-1042 


210 


13.4 2.03 6.45 


MCG -03-12-002 


139 


15.0 0.53 5.51 


IRAS 14378- 


-3651 


294 


14.7 0.93 6.40 


IRAS 16164- 


-0746 


117 


13.2 1.20 5.71 


IRAS 194584 


-0944 


438 


14.4 1.90 7.06 


IRAS 20414- 


-1651 


381 


13.8 1.01 6.66 


'^Luminosity distance calculated using Hq = 70kms ^ Mpc ^ and go = 0.5. 


^1-0 S(l) line flux extracted in 


a 0.75" length of a 1.00" slit; errors are given with the line 


ratios in Table 4. 










Table 3. 


Parameters for the H2 PDR Models 


model X 


nn 


7-" a 
max 


^l-OS(l)^ W 




(cin~^) 


(K) 


(Wm-2sr-i) (Wm-^sr-i) 


1 102 


10^ 


10^ 


1.29 X 10-'' 7.06 X 10-^ 


2 10^ 


10^ 


10^ 


5.17 X 10-^ 2.88 X 10-6 


3 10^ 


10^ 


2 X 10^ 


2.00 X 10-6 3.22 X 10-^ 


4 10^ 


10^ 


10^ 


7.41 X 10-8 4.04 X 10-6 


5 10^ 


10^ 


10^ 


8.82 X 10-^ 1.80 X 10-5 



^For model 1 T = 100 K at all cloud depths. For the other models, the temperature varies as 
given in Eq. 1. 

^The 1-0 S(l) intensity is that predicted by the model. Intensities for other lines can be found 
by using the model line ratios given in Table 4. 

'^The total H2 intensity summed over all the H2 lines is that predicted by the model. 
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Table 4. Relative Line Fluxes 



line and wavelength (/J.m) 



object/model 


iCT* 


1-0S(2) 


He I 
2 0687 


3-2S{5) 


2-lS{3) 
2 0736 


1-0S{1) 
2 1218 


2-lS(2) 
2 1542 


Br7 
2 1661 


3-2S(3) 
2 2014 


l-OS(O) 


2-lS(l) 
2 2477 


MCG -03-04-014 


0.0.53 


0.3.54 


0.316 


<().16() 


0.128 


1.000 


<().16() 


1.254 


<().1B() 


0.298 


0.194 


IRAS 01364-1042 


0.016 


0.365 


0.282 


<0.048 


0.139 


1.000 


0.066 


0.600 


0.062 


0.262 


0.140 


MCG -03-12-002 


0.042 


0.335 


0.221 


<0.126 


0.173 


1.000 


0.078 


0.667 


0.088 


0.218 


0.113 


IRAS 14378-3661 


0.034 


0.325 


0.459 


<0.103 


0.144 


1.000 


<0.103 


1.163 


<0.103 


0.268 




IRAS 16164 - 0746 


0.022 




0.665 


<0.066 


0.169 


1.000 


0.078 


1.331 


0.079 


0.321 


0.138 


IRAS 19468+0944 


0.052 


0.317 


0.386 


<0.167 


0.166 


1.000 


<0.157 


1.038 








IRAS 20414- 1661 


0.061 


0.328 


0.348 


<0.183 


0.113 


1.000 


<0.183 


0.797 


<0.183 


0.268 


0.123 


mean H2 ratios 




0.337 




<0.120 


0.146 


1.000 


0.074 




0.076 


0.271 


0.142 


std. dev. of ratios 




0.018 






0.022 




0.007 




0.013 


0.036 


0.031 


thermal: 1000 K 




0.27 




0.00 


0.00 


1.00 


0.00 




0.00 


0.27 


0.01 


2000 K 




0.37 




0.00 


0.08 


1.00 


0.03 




0.01 


0.21 


0.08 


PDR'' : model 1 




0.49 




0.02 


0.22 


1.00 


0.26 




0.10 


0.64 


0.63 


model 2 




0.32 




0.04 


0.16 


1.00 


0.08 




0.07 


0.33 


0.17 


model 3 




0.35 




0.00 


0.06 


1.00 


0.02 




0.00 


0.22 


0.06 


model 4 




0.31 




0.05 


0.16 


1.00 


0.08 




0.07 


0.32 


0.17 


model 6 




0.29 




0.02 


0.10 


1.00 


0.04 




0.04 


0.33 


0.11 



^The la error applies to all line fluxes given, including 1-0 S(l). It is measured directly from the spectrum as the riTLS of the residual after 
subtraction of both the stellar continuum and emission lines. 

•^For details of the PDR models, see Table 3. 



Note. - 
detected. 



Lines with no ratio given were not within the wavelength covered; 3ct upper limits are given for lines within this range that were not 



Table 5. Derived Data 



object 



log 



log 



Lq 



log 



^Lyc 



b 



QLyc 



100 X 



Qlyc 



MCG -03-04-014 


11, 


,62 


8, 


.83 


53 


.82 


0. 


,33 ± 


0, 


,06 


1, 


,92 


± 


0, 


,06 


IRAS 01364-1042 


11, 


.79 


9, 


.15 


53, 


.68 


0. 


,75 ± 


0, 


,01 


3, 


,58 


± 


0, 


,01 


MCG -03-12-002 


11, 


.48 


8, 


.15 


52, 


.78 


0, 


,86 ± 


0, 


,05 


2, 


,56 


± 


0, 


,10 


IRAS 14378-3651 


12, 


.09 


8, 


.93 


53, 


.91 


0. 


,41 ± 


0, 


,02 


3, 


,03 


± 


0, 


,02 


IRAS 16164-0746 


11, 


.52 


8, 


.72 


53, 


.29 


0. 


,41 ± 


0, 


,01 


3. 


,80 


± 


0, 


,01 


IRAS 19458+0944 


12, 


.27 


9, 


.39 


54, 


.53 


0. 


,48 ± 


0, 


.06 


2. 


,83 


± 


0, 


.06 


IRAS 20414-1651 


12, 


.18 


9, 


.51 


54, 


.02 


0. 


,46 ± 


0, 


.14 


3. 


,33 


± 


0, 


.12 



^Lboi is calculated from the IRAS fluxes for the range 8-1000 //m, using the formula given 
in Tokunaga (2000). 

^Qh2 is calculated from a fluorescently excited H2 hne for PDR model 2. Quoted errors 
include uncertainties in both rates. 
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Table 6. Star Formation Models 



age(Myr) log log ^ ^ 100 x ^ log ^ 



Instantaneous, cluster 


mass 10^ 


Mq, Salpeter IMF, Mupper = 120 M( 


1 


5.48 


51.95 


0.63 


19.02 


8.27 


5 


5.06 


50.78 


5.12 


1.51 


7.89 


10 


5.70 


49.24 


61.85 


0.28 


7.48 


Instantaneous, cluster 


mass 10^ 


M0, Salpeter IMF, M^pp^r = 30 Mq 


1 


4.66 


50.91 


2.25 


9.41 


7.67 


5 


4.63 


50.80 


5.57 


1.05 


7.80 


10 


5.76 


49.31 


61.88 


0.28 


7.55 


Continuous, 


SFR^l 


M0yr-\ 


Salpeter 


IMF, Mupper = 


120 M© 


1 


6.48 


52.95 


0.61 


20.54 


9.25 


3 


6.89 


53.37 


0.80 


16.36 


9.76 


5 


7.03 


53.42 


1.21 


15.10 


9.90 


10 


7.55 


53.43 


1.34 


14.80 


10.03 


100 


8.06 


53.43 


2.02 


14.81 


10.20 


Continuous, 


SFR= 1 




Salpeter 


IMF, Mapper = 


30 Mq 


1 


5.66 


51.91 


2.28 


9.42 


8.66 


3 


6.15 


52.40 


2.89 


8.50 


9.18 


5 


6.37 


52.60 


3.81 


6.00 


9.44 


10 


7.48 


52.66 


5.36 


5.36 


9.74 


100 


8.09 


52.67 


8.66 


5.22 


10.07 
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Fig. 1. — Acquisition images of the targets taken through a Kg filter (with the excep- 
tion of IRAS 01364— 1042 and IRAS 20414— 1651, which were taken through a narrow-band 
2.19/im filter due to the presence of a bright star in the ISAAC field of view). Each field 
is 20"x20", and is drawn with logarithmic scaling. The parallel lines indicate the position 
angle and width of the slit. Due to cosmetic effects, only a narrow strip of the image for 
IRAS 20414-1651 is shown. 
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Fig. 2. — Spectra of each object, plotted at rest wavelength and normalised to their mean 
values. The thick overlaid line is the best fitting continuum constructed from stellar tem- 
plates (see text for details). The vertical marks indicate the positions of the emission lines 
recorded in Table 4. Subtraction of the continuum is necessary in order to detect the faint 
lines such as 3-2 S(3) at 2.2014 /im, which can be seen clearly in IRAS 01364— 1042. 
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Fig. 3. — Excitation diagrams for the mean line ratios, and the 5 PDR models discussed 
in the text. Errors denote the standard deviation in the values used to calculate the mean; 
arrows denote mean of the 3a upper limits. The two dotted lines indicate where the points 
for the best fittinff mirelv thermal single temnerature models would lie: the fits are to all the 
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Fig. 4. — Excitation diagrams for the hot hydrogen molecules in each object. Arrows denote 
3(7 upper limits, derived from the residual spectrum after subtraction of the stellar continuum 
and line emission. The dotted lines indicate where the points for the best fitting purely 
thermal single temperature model would lie. The overlaid grey points show the populations 
for our PDR model 2. 
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Fig. 5. — Upper: Normalised spatial distributions of the continuum flux density (solid line), 
the Br7 flux (dotted hne), and the 1-0 S(l) flux (dashed hne). A single normalisation was 
applied to the Br7 and 1-0 S(l) fluxes. Lower: 1-0S(1)/Br7 ratio (arrows denote 3a limits). 
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Fig. 6. — Position-velocity diagrams. The thick grey hne is the Br7, while the thin black 
line is for 1-0 S(l). The pale vertical stripe shows the region over which the spectra were 
extracted. 



